Hopping Conduction and Bacteria: Transport in Disordered Reaction-Diffusion 
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We report some basic results regarding transport in disordered reaction-diffusion systems with 
birth [A 2A), death (A —* Q), and binary competition {2 A — > A) processes. We consider a model 
in which the growth process is only allowed to take place in certain areas — "oases" — while the rest 
of space — the "desert" — is hostile to growth. In the limit of low oasis density, transport is mediated 
through rare "hopping" events, necessitating the inclusion of discreteness effects in the model. By 
first considering transport between two oases, we are able to derive an approximate expression for 
the average time taken for a population to traverse a disordered medium. 



Reaction-diffusion (RD) models are often used to 
study population dynamics, animal coat patterns, and a 
great variety of other systems, both biological and non- 
biological [1-4]. The effects of quenched spatial disorder 
in the reaction rates on the behavior of these models 
have been difficult to discern — in some cases, a renor- 
malization group analysis yields runaway flows [5] — but 
some progress has been made [4, 6-9]. In particular. 
Nelson and coworkers have looked at the effects of con- 
vection and quenched spatial disorder on the evolution 
of a population density described by a generalization of 
the Fisher/KPP equation given by 



dc{x, t) 
dt 



D\7^c{x, t) + U{x)c{x, t) - qc(x, tf , (1) 



where c{x, t) represents the population density, Z? is a 
spatially homogenous diffusion constant, U{x) is a spa- 
tially inhomogeneous growth term, and q — Mq'^ is a com- 
petition term (6 is a competition rate and £a is the mi- 
croscopic length scale at which two particles will compete 
with one another) [10]. One simple form of inhomogene- 
ity considered in these works is a "square well" potential 
U{x) which consists of a uniform space with negative 
growth rate — termed the "desert" — in which a single re- 
gion of positive growth rate — an "oasis" — is placed. This 
model has proven to be applicable to experiments with 
bacteria populations in adverse environments [11]. 

Other important studies have focused on the effects 
of spatial inhomogeneities in both reaction and diffu- 
sion rates on the speed and roughness of front solutions 
to reaction-diffusion equations [12, 13]. These works 
consider a parameter regime — which we dub fertile — in 
which a front solution exists in the absence of disorder 
(for (1), this regime corresponds to U{x) = U > 0). The 
nature of transport in the opposite parameter regime in 
which the average growth rate is negative — which we dub 
the hostile regime — has not been studied as much. In 
this regime, there may be localized stationary solutions 
centered around regions of positive growth rate — oases 
[10] — as well as noise-induced front solutions for suffi- 
ciently large variations in growth rate [14]. 

In this letter, we report some results on the nature of 



transport in a hostile disordered system; specifically, we 
consider the case of a desert into which identical oases 
are placed randomly at low density. Because transport 
between oases involves the movement of a low popula- 
tion density, it is natural to assume that discreteness 
effects — fluctuations about the mean field theory — might 
come into play. For this reason, we consider the stochas- 
tic dynamics underlying the differential equation (1) in 
which discrete particles are allowed to diffuse, reproduce 
{A -> 2A), die {A^ 0), and compete {2 A A). Using 
a novel method, we find the probability density function 
for the time at which a particle from a populated oa- 
sis first reaches an unpopulated oasis — the first passage 
time; we then show how, for sufficiently low oasis density, 
this result can be used along with the theory of hopping 
conduction to estimate the mean transit time across the 
system. 

We start by studying transport between two oases, 
one of which — we will call this the first oasis — is initially 
populated. We define the infection time to be the time 
elapsed before the second oasis reaches the population 
level at which the first oasis started. This time can be 
broken into two parts: the time Jtransit that it takes the 
population to reach the second oasis and the time Tgrowth 
it takes for the population at the second oasis to grow to 
the specified level. We assume very fertile oases; that is, 
we assume that the particles which first reach the second 
oasis will reproduce immediately, ignoring the possibility 
that they may return to the desert and die. (This can 
be accomplished by having oases with high growth rates 
or by "seeding" the oases with a second species B and 
including a new very fast reaction A + B 2A.) This 
makes Ttransit identical to the first passage time (FPT) 
of the process — that is, the time it takes for the first 
particle to reach the second oasis. We will calculate the 
probability distribution of this time — the first passage 
time probability distribution function (FPT PDF). 

Consider two oases of radius a in c? dimensions whose 
centers are separated by a distance R. The oases each 
have growth rate y, and the desert between them has 
death rate z. The competition rate is Mq'^, where Iq is a 
microscopic length scale. If R is sufficiently small, we ex- 
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pect that the FPT should be hnear in R: a "wavefront" of 
particles radiates out from the first oasis due to diffusion, 
its amplitude decreasing due to the death term, eventu- 
ally reaching a level small enough that mean field theory 
is invalid. The distance i?iin at which this linear transport 
breaks down, then, should be approximately the same as 
the distance at which mean field theory breaks down. 
This distance can be obtained by solving the mean field 
differential equation (1) and looking at the solution for 
large times; the distance where the population density 
falls off to ~ io"' is the distance we are looking for. In 
one dimension, c{x) ^ cq e~\f^^^ for large x. The con- 
stant Co is of order y/{MQ^) for small values of y, and 
thus an estimate for is \jD/z \n{y/b). In higher di- 
mensions, the relevant length scale is smaller, but we will 
use the above expression as a rough estimate. 

Because the competition process 2A^ A\s not a sig- 
nificant mechanism of particle destruction for distances 
beyond i?iin, the particles' interactions with each other 
can be ignored. With this simplification in mind, a model 
for the behavior of the system at large oasis separa- 
tion can be formulated by replacing the first oasis with 
desert and a point source (0 A) which creates N non- 
interacting particles per unit time. If the creation rate is 
chosen so that the average flux of particles at _Riin is the 
same as for the model with competition, the FPT statis- 
tics should be nearly identical for large oasis separations. 
We will refer to this simplified model as the linear model 
with a source, and to the full model with competition as 
the nonlinear model. 

Since the particles are non-interacting in the linear 
model with a source, the full FPT PDF, which we will 
denote /jv(-R, i), can be written in terms of the one- 
particle FPT PDF fi{R,t). We will do this as follows: 
let the source be at the origin, and define S{R, t) = 
1 - J^^dt' fi{R,t') = 1 - Phit{R,t) to be the probabil- 
ity that a particular particle has not reached the oasis 
located a distance R away by time t. The probability 
Pnone{R,t) that this sitc has never been visited by any 
particle by time t is just a product of the probabilities 
S{R, t) that each of the particles has never visited it: 



,{R,t) 



n 

T=0,At,. 



[S{R,r)] 



N 



(2) 



Taking the natural logarithm, letting At ^ 
with N/At = g fixed, and re-exponentiating gives 

Pnone{R,t) = cxp (jj Jq dt' liiS{R,t')j . Sincc we are in- 
terested in sites sufficiently far from the origin {R » 
y/D/z) that a given particle has low probability of ever 
reaching, we can approximate In 5 = ln(l — Phit) by 
— fhit, which leads to the following useful expression for 



{R, t) ~ exp 



1 [ dt'{t- 
Jo 



The FPT PDF fN{R, t) is related to Pnone by /jv(ii, t) = 

— dtPnone{Rj t) . 

We will apply these ideas first to a ID system on a 

lattice without convection. On a lattice, the diffusion 
constant D is replaced by a hopping rate w (the rate 
for hopping to a particular side is w/2). We replace the 
first oasis centered at the origin with a source and ask 
for the FPT PDF to the i/-th lattice point for the lin- 
ear model. By using the proper one-particle FPT PDF 
fiii^A) = |z/|e-("'+^)*/j.(u;t)/t [15] along with (3), we 
can find fN{i^, t) and all of its moments. In order to com- 
pare these values to simulations of the nonlinear model, 
it is essential to fix g to an appropriate value. To do 
this, we match the long-distance t ^ co solutions for the 
mean particle concentration in the linear model with a 
source and the nonlinear model. The linear model has a 
t ^ oo solution of Coojinlj^) = ge~^\'^^/'wsmh{f ), where 
/ = cosh~^(l + z/w), while the nonlinear model decays 
like 

Coc.nlin 

(ly) ~ coe ■'I'"!. We have determined co and 
thus g numerically. 

For the d = 1 lattice case, fN{i^, t) is given by a com- 
plicated expression, but it does have a simple asymptotic 
behavior, decaying exponentially like e^^^''^* as ^ ^ oo, 
with = gie"-^!'^!. The j-th moment of /jv(i', — 

which we denote {T^{v)) — approaches a simple asymp- 
totic limit as v ^ oo [16]: 



(ID lattice) 



(4) 



t')MR,t') . (3) 



The first moment of fN{i^,t) — the mean first passage 
time thus depends exponentially on the separation of 
the oases in the limit of large oasis separation. 

To test the predictions made using the linear model 
with a source, we wrote a kinetic Monte Carlo program 
to simulate the full discrete stochastic process in one di- 
mension. The agreement between the Monte Carlo re- 
sults and the linear theory with a source is excellent. 
The linear theory correctly predicts the mean FPT for 
distances sufficiently far from the oasis, as can be seen in 
Fig. 1. An even more stringent test of the power of the 
linear theory is a comparison of its prediction of the full 
FPT PDF with simulation; this is shown in Fig. 2. Even 
with only 5000 nms, the simulation data begins to fill 
out the shape of the FPT PDF predicted by the linear 
theory with a source. 

With g specified, the only remaining source of ambi- 
guity in matching up the predictions of the linear theory 
with a source to the full nonlinear simulation results is 
the choice of simulation initial conditions. This ambi- 
guity stems from the fact that the early time dynamics 
of the linear theory and the nonlinear theory are differ- 
ent; they take different times to "grow up" to the point 
where their mean particle fluxes into the desert are equal. 
This disparity becomes less important as transit to sites 
further away is considered. For the simulations we per- 
formed, we chose to place y/2b = 125 particles (half the 
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Monte Carlo Data i 



of this distribution are given in the large R hmit by [16] 



n ( latt ice point ) 



FIG. 1: Comparison of theoretical mean FPT with Monte 
Carlo data. The error bars represent a 95% confidence inter- 
val. The parameters used were w = 1.0, y = .25, z = .1, and 
b = .001. Time is measured in units of 1/w. 
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FIG. 2: Comparison of theory with Monte Carlo data. Each 
box represents the probabihty P{ti) (from the linear theory) 
that a particular run hits site n = 25 for the first time in the 
i-th time bin. The points with error bars are from 5000 Monte 
Carlo simulations of the full nonlinear stochastic model. Time 
is measured as in Fig. 1, and the bins have a width of 25 units 
of time. 



carrying capacity) at the center site of the oasis of 5 total 
sites at t — 0, with no particles at any other site. This 
choice led to first passage times in agreement with the 
linear theory for the parameter sets considered. 

The continuum case in one dimension is actually sim- 
pler to deal with than the lattice case previously de- 
scribed, even when a uniform convection velocity v is 
included. (Physically, v is an external influence on the 
system, like the motion of a liquid in which bacteria 
live.) The FPT PDF /Ar(i?,i) again has a very com- 
plicated form, but its behavior in certain limits is eas- 
ily described. For any R, the function decays exponen- 
tially as e"^^"'*-'* as t ^ oo, with 7(-R) = g{v) e""^-'*-) and 
a{R) = v^l + /4Dz ^/z/DR- vR/2D. The moments 



(T^(i?))=j! 



(ID continuum) 



(5) 



For higher dimensions (d > 1), the lattice problem be- 
comes quite complicated, while the continuum problem 
is still analytically tractable. Instead of calculating the 
single particle FPT PDF fi{R,t) to a point, one instead 
finds fi{R,a,t) to a hypersphere of radius a, given an 
initial condition in which the particle is located at posi- 
tion R a distance \R\ = R > a from the center of the 
hypersphere [15]. As in one dimension, this function is 
then used along with (3) to determine fN{R,a,t). For 
zero convection velocity, we have obtained the following 
result for the moments of the FPT PDF as i? — > oo [16]: 



{T^{R,a)) =j! 



(6) 



where is the /i-th order modified Bessel function of 
the second kind. As in the d — 1 continuum case, it 
is possible to include the effects of convection. For a 
small uniform convection velocity \v\ <C \JDz^ the 
j-th moment is multiplied to leading order by a factor 
\g(v = ^)e^'^^'^^ I g[vf\ . Previous studies concentrat- 
ing on (1) suggest that g(v) should decrease as \v\ in- 
creases [10]. 

We have shown that a linear model with a source can 
capture the long-distance behavior of the FPT PDF for 
a nonlinear model with two oases. We wish to apply 
our results to a system with many oases. Consider a 
very large d-dimensional (d > 1) continuum system in 
which identical oases of radius a and birth rate y are 
centered (with overlaps allowed) around randomly placed 
points with number density ti in a desert of death rate 
z. (Grassberger has studied the related problem of ran- 
domly placed traps in a neutral background [17], and 
Redner has studied lattice systems with mixed oases and 
traps [18].) We are interested in the low oasis density 
regime in which the average distance between oases is 
much larger than i?iin- One oasis is populated at i = 0; 
as time goes on, the population will spread , and eventu- 
ally an oasis located at position L will be reached. We 
wish to find the infection time Tinfcction — the time for this 
process to occur. 

Because of the exponential dependence of Kii/2-i{KR) 
on R for large R in (6), there is a wide distribution of 
mean transit times between oases in the system. The sit- 
uation is mathematically analogous to (though physically 
quite distinct from) hopping conduction in doped semi- 
conductors: the oases are the analogs of impurity sites, 
and the mean transit time plays the role of the resistance 
between sites, which is roughly equal to e"^/Go, where 
Go is a constant. In the hopping conduction problem, 
the resistance of the system is dominated by the pairs of 
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impurity sites with the largest separations in the cluster 
(the critical subnetwork) which carries the bulk of the 
current [19, 20]. This cluster has correlation length Lq, 
and above this length scale the system can be regarded 
as roughly homogenous. The resistivity p is then given 
by p ~ Lo'^~^e"^'°»''/Go, where iimax is the largest dis- 
tance between impiirity sites in the critical subnetwork, 
which is roughly equal to the critical percolation radius 
Rc obtained from continuum percolation theory [19]. It 
should be noted that there is some debate as to whether 
Lq is the proper distance to use in d = 3 for the typi- 
cal separation of large resistances; one competing theory 
identifies an additional length scale I which enters into 
the resistivity along with Lq [21]. 

In our problem, the largest oasis separations should 
dominate the transit time statistics across a sample of 
size Lq. This suggests a coarse-graining that can be used 
to determine the infection time: we replace each block of 
size Lo'^ with a node, and assign to each node a random 
transit time picked from /Ar(i?max, a, i), the FPT PDF 
for crossing the largest oasis separation in each node. 
This is a rough approximation; the important point is 
that the mean time to cross each block is largely deter- 
mined by the largest oasis separations. We ignore block- 
to-block variations in the size of the largest oasis sepa- 
ration, since these become small at the length scale Lq 
[22] . The infection time is now given by the path with the 
shortest total transit time that goes from the block/node 
containing the starting oasis to the block/node contain- 
ing the target oasis; we have turned our problem into a 
first passage percolation (FPP) problem [23] . It has been 
shown that, as the distance between nodes goes to infin- 
ity, Tinfection/?7i, where m is the distance between starting 
and target nodes, approaches a constant /x, convention- 
ally called the time constant [23]. An upper limit for jj, 
is given by the mean time to cross one node. In the ap- 
proximation we have made, this is (T(-Rmax. o))- Thus, 
an estimate for the mean infection time to a distant tar- 
get site at L divided by |X| is: 

(^Infection) {T {Rmax t d)) /_-, 

\L\ - To ' 

where i?max is given by continuum percolation theory as 
-Rmax = [Bc{d)/{nVd)Y^'^, where Va is the volume of a d- 
dimensional unit hypersphere and Bc{d) is a dimension- 
less number known as the bonding criterion. The value 
of (T(i?max5 a)) is given approximately by (6) for large 
R. Lq is roughly equal to [-^/zJD R^^xY , where v 
is a critical exponent (~ 4/3 in = 2, .88 in li = 3) [19]. 

We have presented both analytical and numerical re- 
sults of investigations of transport in disordered reaction- 
diffusion systems. Taking discreteness effects into ac- 
count, we have shown that the first passage time between 
two oases separated by a distance i? in a desert is a very 
broadly distributed quantity whose mean increases expo- 
nentially with R for large R (see Eq. 6). We have used 



an analogy with hopping conduction to argue that the 
largest oasis separations dominate transit times up to a 
length scale Lq, and we have employed a mapping to a 
first passage percolation system to then arrive at an es- 
timate (7) for (Tinfection)/!-^^!, the mean infection time 
to a site located at L divided by \L\. We are currently 
working on determining the effects of a uniform convec- 
tion velocity on our model as well as the nature of the 
advancing front for the case where half of space, rather 
than one oasis, is initially populated. 
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